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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  26^1 

A  VECTOR  STUDY  OF  LINEARIZED  SUPERSONIC  FLOW 

APPLICATIONS  TO  NONPLANAR  PROBLEMS 


By  John  C .  Martin 


SUMMARY 


A  vector  study  of  the  partial-differential  equation  of  steady 
linearized  supersonic  flow  is  presented.  General  expressions,  which 
relate  the  velocity  potential  in  the  stream  to  the  conditions  on  the 
disturbing  surfaces,  are  derived.  In  connection  with  these  general 
expressions  the  concept  of  the  finite  part  of  an  integral  is  discussed. 

A  discussion  of  problems  dealing  with  planar  bodies  is  given  and 
the  conditions  for  the  solution  to  be  unique  are  investigated. 

Problems  concerning  nonplanar  systems  are  investigated,  and  methods 
are  derived  for  the  solution  of  some  simple  nonplanar  bodies.  The  sur¬ 
face  pressure  distribution  and  the  damping  in  roll  are  found  for  rolling 
tails  consisting  of  four,  six,  and  eight  rectangular  fins  for  the  Mach 
number  range  where  the  region  of  interference  between  adjacent  fins 
does  not  affect  the  fin  tips. 


INTRODUCTION 


In  the  presentation  of  the  theory  of  the  flow  of  an  idealized 
incompressible  fluid,  vector  methods  can  be  used  to  reduce  greatly  the 
mathematical  manipulations  involved.  The  study  of  steady  linearized 
supersonic  flow  may  also  be  aided  by  the  use  of  vector’ methods .  Two 
types  of  approaches,  however,  can  be  used.  Perhaps  the  more  obvious  is 
to  make  use  of  common  vector  methods  as  was  done  in  reference  1.  The 
other  vector  method,  which  was  introduced  by  Robinson  in  reference  2 
and  is  used  in  this  paper,  appears  to  be  more  suited  to  the  study  of 
the  linearized  partial-differential  equation  of  steady  supersonic  flow. 
This  method  allows  a  derivation  of  a  hyperbolic  scalar  potential  and  a 
hyperbolic  vector  potential  along  lines  analogous  to  the  derivation 
sometimes  used  (reference  3*  ch.  VIII)  in  dealing  with  common  scalar 
and  vector  potentials. 
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The  present  paper  presents  a  vector  derivation  of  many  general  4 

results  which  have  been  found  by  various  methods  and  are  given  in  the 
published  literature  on  the  linearized  partial-differential  equation 
of  supersonic  flow  and  also  presents  some  results  which  are  not  found  v> 

in  the  literature.  The  general  results  of  Hadamard  (reference  4, 
p.  207)>  Puckett  (reference  5),  and  Heaslet  and  Lomax  (reference  6) 
are  fpund  as  special  cases  of  a  general  expression  for  a  scalar  poten¬ 
tial,  and  the  results  found  by  Robinson  (reference  2)  are  obtained  by 
the  use  of  a  vector  potential.  The  derivation  of  the  scalar  potential 
doubtlessly  helps  to  clarify  the  concept  of  the  finite  part  of  an 
integral . 

A  discussion  of  problems  dealing  with  planar  bodies  immersed  in  a 
supersonic  flow  is  given,  and  the  conditions  necessary  for  the  solution 
to  be  unique  are  investigated. 

Problems  dealing  with  nonplanar  systems  are  also  discussed,  and 
methods  are  derived  for  the  solution  of  some  simple  problems  dealing 
with  nonplanar  bodies.  The  surface  pressure  distribution,  the  span- 
wise  loading,  and  the  damping  in  roll  are  found  for  rolling  tails  con¬ 
sisting  of  four,  six,  and  eight  rectangular  fins  for  the  Mach  number 
range  where  the  region  of  interference  between  adjacent  fins  does  not 
effect  the  fin  tips. 

SYMBOLS  ,  » 


A 

A 

a 

C1,C2,C3 

C 

E,F 

f 

G 

H 

b/2 


hyperbolic  vector  potential 
aspect  ratio  of  tail  fin 
positive  constant 
arbitrary  constants 

chord 

arbitrary  vector  functions 

scalar  function  defined  by  equation  (19) 

vector  function  associated  with  vector  function  F 

vorticity  vector 

span  of  tail  fin 


t 
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-unit  vectors  in  x-,  y- ,  and  z-directions,  respectively 
M  Mach  number 

n  unit  vector  normal  to  element  of  area  da 

nh  =  -ip2V1  +  jv2  +  kv3 


r^’  =  -ip2V^ '  +  jv2'  +  kv3' 
Bj*  =  +  jv2*  +  kv3* 


q' 


pres sure -difference  coefficient 
rate  of  roll 

function  used  in  equation  of  surface  of  discontinuity 

part  of  velocity  vector  which  is  made  up  of  hyperbolic 
curl  of  vector  potential 

total  perturbation  velocity 


R  =  /(x  -  |)2  -  p2(y  -  ti)2  -  P2(z  -  02 
R1  small  constant 


r  =  \J(x  -  t)2 
S 

S0 

SijS2,s3,s^, 

s5,s6,T,T' 

V 

v0-’ vl,v2 


+  p\y  -  q)2  +  P^Cz  -  02 

area  of  tail  fin 
surface  of  discontinuity 

surfaces  of  integration 

free -stream  velocity 
volumes  of  integration 


W  -  ^ 

x,y,  z  Cartesian  coordinates  (x-axis  parallel  to  free-stream 

direction) 
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=  sfiF 


r 

€ 

6 

e,p 

\ 


V1,V2,V3 
V1  ,V2  >v3 


v  *  v  *  v  * 
V1  >v2  ’v  3 


spanwise  circulation  ^2 
small  positive  quantity 


rT‘E-  ,<  „  V  rT*E- 

/  0x  dx  =  p  /  ACp  dx 

L.E.  d  JL.-E. 


Cartesian  coordinates  (g-axis  parallel  to  free-stream 
direction) 

polar  coordinates 

scalar  functions 

given  volume 

direction  cosines  of  outward  normal  to  element  of  area  da 

direction  cosines  of  normal  (directed  away  from  point 
(x,y,z))  to  surface  SQ 

direction  cosines  of  normal  to  element  of  area  da  used 
in  equation  (45) 


a  slope  of  deflected  area 

t  area  of  integration 

0,00^1^2^0  scalar  potentials 


rolling-moment  coefficient  per  fin 

(Rolling  moment  per  fin\ 

) 


Cl  = 

Lv 


6c 


p(V2) 


p— >0 


indicates  integration  over  closed  line  or  surface 


denotes  finite  part  of  integral 


r 
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THEORY 


This  paper  deals  with  the  linearized  partial- differential  equation 
of  steady  supersonic  flow.  This  equation  is  given  by- 

-P20  +  0  +  Of  =  0  fi) 

rxx  ryy  rzz 

The  potential  is  assumed  to  he  continuous  in  the  stream  direction,  and 
the  potential  is  assumed  to  he  always  finite.  Assuming  the  potential 
to  he  finite  and  continuous  in  the  stream  direction  has  the  effect  of 
requiring  the  aerodynamic  lift  and  moment  (calculated  hy  use  of  the 
linearized  pressure)  of  finite  bodies  to  he  finite  since  the  linearized 
pressure  is  related  to  the  derivative  of  the  potential  in  the  stream 
direction.  The  expression  linearized  pressure"  refers  to  the  pressure 
obtained  hy  neglecting  all  powers  of  the  perturbation- velocity  compo¬ 
nents  above  the  first. 


Vector  Operators  and  Identities 

Certain  operators,  which  are  closely  associated  with  the  linearized 
hyperbolic  partial-differential  equation  of  supersonic  flow  (the  two- 
dimensional  wave  equation) ,  are  added  to  the  vector  operators  commonly 
used.  The  basic  operators  have  been  used  previously  in  references  2 
and  7- 

The  gradient  operator  is  defined  by 


v  =  1  !e  +  +  k  Jb 

The  analogous  hyperbolic  gradient  operator  defined  by  Robinson  in 
reference  2  may  be  expressed  as 


Vh  -  -ip2  J-  +  j  +  k  |_ 

The  hyperbolic  divergence  of  an  arbitrary  vector  E  is  given  by 

Vh  •  E 

Similarly,  the  hyperbolic  curl  of  the  vector  E  is  given  by 
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The  divergence  of  the  gradient  operator  is  sometimes  denoted  by 


The  analogous  divergence  of  the  hyperbolic  gradient  operator  is  denoted 

by 


Vh 


The  following  identities  are  needed.  Let  E  be  a  vector  and  i  and  A 


be  scalar  functions  of  x,  y,  and  z.  Then, 

Vi  •  VA  =  VA  ■  Vi  (2a) 

V  •  \|/E  =  iV  •  E  +  E  •  Vi  (2b) 

V  X  (V  X  E)  =  V(V  •  E)  -  V%  (2c) 

V  •  (V  X  E)  =  0  '  (2d) 

Vhi  •  VA  =  VhA  •  Vi  (2e) 

Vh  •  i|rE  =  iVh  •  E  +  E  •  Vhi  (2f) 


V  X  (Vh  X  E)  =Vh(v  •  E)  -  V2hE 


(2g) 


Vh  X  (V  X  E)  =  V (Vh  •  E)  -  V2hE 
Vh  •  (vh  X  E)  =  0 

These  identities  can  he  proved  hy  direct  expansion. 
The  divergence  theorem  may  he  expressed  as 


n  da 


V  •  E  dv 


(2h) 

(2i) 


(3) 


where  n  is  the  normal  unit  vector  to  the  element  of  area  da.  The 
vector  n  is  expressed  mathematically  as 
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where  V and  are  the  direction  cosines  of  the  outward  drawn 

normal  to  the  element  of  area  da. 


A  theorem  more  general  than  the  divergence  theorem  is  given  by 
(this  theorem  follows  from  the  results  of  reference  8,  p.  87) 


/(°i 


vlEx 


CgVgE. 


C3V3Ez 


)da  -/(c. 


3ex 

sr 


+  c 


cSEy 

2  5^ 


+  C. 


dEz 

5z“ 


•W 


/ 


where_jthe  subscripts  x ,  y,  and  z  refer  to  components  of  the  vec¬ 
tor  E ,  and  C2,  and  are  arbitrary  constants.  Note  that  if 

=  1  the  preceding  equation  reduces  to  equation  (3).  If 


C 


1 


-  -P 


2 


C 


2 


=  1 


the  preceding  equation  reduces  to 


H- 


2 

P  VnE 
lx 


V0E 
2  y 


V0E 
3  z> 


da 


BEX  SEy 

- -  +  -  4- 

ox  ay 


or 


n^  da 


where 

nh  =  +  jv2  +  kv^ 

If  the  divergence  theorem  as  expressed  hy  equation  (3)  is  applied 
to  a  volume  throughout  which 

V  •  E  =  0  * 

then  the  surface  integral  over  the  bounding  surface  is 
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provided  that  no  surfaces  exist_inside  the  volume  of  integration  across 
which  the  normal  component  of  E  is  discontinuous.  Similarly,  if 
equation  (4)  is  applied  to  a  volume  throughout  which 

Vh  •  E  =  0 

then  the  surface  integral  over  the  hounding  surface  is 


nh  da  =  0 


provided  that  there  are  no  surfaces  inside  the  volume  of  integration 
across  which  E  •  n^  is  discontinuous.  It  is  interesting  to  note, 
however,  that  surfaces  exist  inside  the  volume  of  integration  across 
which  E  *  n  can  he  discontinuous  while  at  the  same  time  E  • 
remains  continuous.  It  follows  that  for  such  a  surface  n  and  n^ 
must  satisfy  the  relation 


n  •  nh  =  0  (5) 

Let  Q(x,y,z)  =0  be  the  equation  of  such  a  surface. 

Then , 

n  =  -  -  h  : - -  VQ 

\j%?  +  Qy2  +  Qz2 

and 


n 


h 


where  the  subscripts  indicate  differentation.  Substituting  the 
preceding  expressions  for  n  and  into  equation  (5)  yields 


* 


(6) 


Any  solution  of  equation  (6)  set  equal  to  zero  is  the  equation_of  a 
surface  across  which  V  •  E  may  he  discontinuous  while  Vh  •  E  remains 
continuous.  The  fact  that  the  Mach  cone  from  any  arbitrary  point  satis¬ 
fies  equation  (6)  can  he  easily  verified.  The  equation  of  the  envelope 
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of  the  Mach  cones  from  an  arbitrary  line  also  satisfies  equation  (6) 
(reference  9,  p.  106). 


Finite  Part  of  Integrals  Which  Arise 
in  Steady  Supersonic  Flow 

In  the  following  sections  use  is  made  of  the  concept  of  the  finite 
part  of  an  infinite  integral.  This  concept  was  introduced  by  Hadamard 
(reference  4)  and  has  been  used  by  a  number  of  other  investigators. 

The  finite  part  is,  however,  sometimes  confusing.  This  section  was 
therefore  included  in  an  attempt  to  give  a  realistic  picture  of  the 
finite-part  concept #and  also  to  present  the  first  steps  of  the  deriva¬ 
tion  of  the  scalar  and  vector  potentials. 

The  concept  of  the  finite  part  of  double  integrals  as  defined  by 
Hadamard  and  used  in  this  paper  is  different  from  the  concept  of  the 
finite  part  of  double  integrals  as  defined  in  reference  10.  The  essen¬ 
tial  difference  between  these  two  definitions  lies  in  the  manner  in 
which  the  singular  points  along  the  Mach  cone  are  treated.. 

In  reference  3,  page  183,  a  vector  function  is  used  in  the  deriva¬ 
tion  of  the  common  scalar  and  vector  potentials.  The  analogous  vector 
function  based  on  equation  (l)  is 

w  =  i  V0  -  fa  i 

where 

R  =  \/(x  -  If  -  p2(y  -  T| ) 2  -  P2(z  -  02 

The  hyperbolic  divergence  of  vector  W  with  respect  to  variables 
q,  and  £  is  given  by 


The  preceding  equation  indicates  that  the  hyperbolic  divergence  of  the 
vector  W  set  equal  to  zero  yields  the  partial-differential  equation 
of  linearized  supersonic  flow.  A  mathematical  derivation  of  W  can  be 
obtained;  however,  for  the  purposes  of  this  paper  such  a  derivation  is 
not  needed. 
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The  result  of  applying  equation  (4)  to  the  vector  W  is 


nh  da  = 


i-V2^  dv 
R 


(7) 


When  0  satisfies  equation  (l)  throughout  the  volume  of  integration, 
the  right  side  of  equation  (7)  is  zero;  thus. 


Z(iv0  -0vi)  .  Shda  = 


0 


(8) 


when 


V2hjZ i  =  0 

Equation  (7)  is  applied  to  a  volume  (denoted  by  Vq)  enclosed  in 
the  forward  Mach  cone  from  the  point  (x,y, z).  This  volume  is  bounded 
by  the  surface  given  by  R  =  R  ’  ,  where  R1  is  a  small  constant,  and  an 
arbitrary  surface  S]_  enclosed  in  the  forward  Mach  cone  from  the  point 
(x,y,z).  A  cross  section  of  the  region  of  integration  is  shown  in  fig¬ 
ure  1.  Note  that  this  region  is  analogous  to  the  region  that  is  some¬ 
times  used  in  calculating  the  potential  function  satisfying  Laplace Ts 
equation  (reference  3,  pp.  151-153)-  For  regions  such  as  the  one  shown 
in  figure  1,  equation  (7)  may  be  written  as 

/T  |  V0  -  07  i)  .  nh  da  +  f  ^  .  nh  da  -  f  i  V2h0  dv 


(9) 

where  T  represents  the  area  of  integration  when  R  =  R ' . 

The  integral  over  the  area  T  may  he  reduced  to 

5^  •  Eh  +  <10> 

where  r  is  given  by 

r  =  fiT-  I)2  +  p\y  -  t^)2  +  p\z  -  02  (11) 


r 


i 
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Since  R1  is  a  constant,  equation  (10)  can  be  written  as 


Equation  (9)  can  now  be  written  as 

•  Hh  +  T^)da  +  ^  (e  *  -  ^l)  -'I4*  -/To  B  ^  4T 

(12) 

If  0  is  required  to  satisfy  the  linearized  partial-differential 
equation  of  steady  supersonic  flow,  then 

V2h0  =  0 

and  equation  (12)  reduces  to 

&f-  ^  +  fi^)da  +/  (l  Vj<  -  *7  i)  .  5h  to  -  0  (13) 

If  R'  is  made  smaller  and  smaller  the  integrand  of  the  integral  over 
the  area  T  in  equation  (13)  remains  finite  except  on  the  small  area 
close  to  the  point  ( x,y,z ).  In  anticipation  of  taking  the  limit  of 
equation  (13)  as  R'  approaches  zero,  the  small  area  close  to  the 
point  (x,y,z)  is  removed  from  the  area  T.  The  area  T  is  divided 
into  two  parts.  One  part  is  the  area  of  T  which  is  downstream  of  the 
surface  given  hy 


I  =  x  -  e 

where  €  is  small  hut  larger  than  R' .  This  area  is  denoted  hy  t. 
The  remaining  part  of  T  (denoted  hy  T')  is  the  area  of  T  which  is 
upstream  of  the  surface 

I  =  x  -  e 

A  cross  section  of  the  region  of  integration  with  T  divided  into  t 
and  T'  is  shown  in  figure  2.  Equation  (13)  can  now  he  expressed  as 
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^  +^4  (^v0- 

/  (s  ^  ^  s)  ■  5h  da  ■ 


n^  da  + 


where  R’  is  smaller  than  e. 

Since  f  is  continuous  and  therefore  its  values  over  t  are 
approximately  constant  for  small  values  of  c  ?  the'  integral  over  the 
area  t  can  he  written  as 


V0  .  n^  da  + 


P20(x,y,z 

RT 


When  the  second  integral  of  expression  ( 15)  is  integrated,  equation  (l4) 
becomes 


V0  .  nh  da  *  .  2*(w)  ♦  i/ >  +  £ 


i  V0  -  0V  .  nh  da  =  0 


If  R'  is  made  to  approach  zero,  equation  (l6)  applies  even  to  the 
limit  where  Rl  is  zero. 

The  limit  of  equation  (l6)  as  R’  approaches  zero  may  be  written 


a,%  V*  -  \  da  +  +  §?/  i  da  ♦ 


/  (iv0  -  0vi)  -  n„  d^  .  0 


The  integrands  of  the  integrals  over  the  areas  t  and  T'  are  always 
finite  and  it  can  he  shown  that  their  first  derivatives  with  respect 
to  R'  approach  zero  as  R'  approaches  zero;  therefore,  the  product 
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of  i/r'  and  these  integrals  either  approaches  zero  in  at  least  the 
order  of  R'  or  approaches  infinity  as  R'  approaches  zero.  Thus  it 
v  follows  that  the  integrals  over  the  areas  t  and  T'  have  no  finite 

terms  remaining  after  the  limit  (R'  — >  0)  has  been  taken.  The  Siam  of 
the  terms  of  equation  (17)  must  he  zero;  thus  the  singularities 
resulting  from  the  integrals  over  the  areas  t  and  T'  must  cancel 
the  singularities  which  arise  from  the  integral  over  the  area  Sq. 

From  the  preceding  considerations  it  follows  that  one  method  of 
evaluating  the  finite  part  of  infinite  integrals  of  the  type  appearing 
in  equation  (17)  is  to  evaluate  the  integral  when  R'  is  small  hut  not 
zero  and  neglect  the  terns  multiplied  by  powers  of  l/R' .  Other  infinite 
integrals  sometimes  arise,  however,  for  which  the  finite  part  cannot  be 
obtained  by  neglecting  powers  of  l/R' .  For  example,  if  equation  (17) 
is  differentiated  with  respect  to  one  of  the  variables  (x,  y,  or  z) 
an  equation  containing  the  velocity  component  is  obtained.  In  some 
cases,  when  the  point  (x,y,z)  lies  on  the  surface  Sq  the  infinite 
terms  are  of  the  order  (in  R')/R'  and  of  the  orders  (l/R')n.  In  these 
cases,  the  finite  part  of  the  infinite  integrals  can  be  obtained  by 
evaluating  the  integrals  when  R'  is  small  and  neglecting  the  terms 
multiplied  by  powers  of  l/R'  and  (in  R')/R'. 

The  process  of  removing  the  infinite  parts  of  an  integral,  however, 
has  been  derived  by  Hadamard  (reference  4,  book  III,  ch.  I).  Hadamard 
used  his  methods  of  evaluating  the  finite  part  of  integrals  in  finding 
solutions  to  certain  hyperbolic  equations  including  the  linearized  equa¬ 
tion  of  steady  supersonic  flow.  Perhaps  a  fact  worth  noting  is  that 
the  integrals  of  equation  (17)  are  double  integrals  and  when  the  methods 
given  by  Hadamard  are  used  the  methods  given  for  multiple  integrals 
should  be  used.  In  the  past,  the  singular  points  (points  on  the  Mach 
cone  where  the  derivative  of  (x  -  |)2  -  p2(y  -  t))2  -  p2(z  -  £)2  with 
respect  to  the  variable  of  integration  is  zero)  have  caused  some  con¬ 
fusion;  as  Hadamard  points  out  (reference  4,  p.  147),  these  singular 
points  must  be  removed  from  the  area  of  integration  before  the  finite 
part  is  taken.  Particular  attention  should  be  given  to  paragraph  92 
of  reference  4  since  the  special  type  of  integrals  discussed  therein 
sometimes  arises  in  dealing  with  planar  problems. 

Robinson  (reference  2)  has  shown  that  when  using  Hadamard' s  methods 
the  order  of  integration  may  be  changed  without  affecting  the  finite 
part  and  that  it  is  permissible  to  differentiate  under  the  integral 
sign  of  a  multiple  integral  without  considering  the  variable  limits 
which  lie  along  the  boundary  where  the  integrand  is  singular,  provided 
that  only  the  finite  part  is  taken.  Both  Hadamard  and  Robinson  have 
shown  that  in  differentiating  an  improper  integral  which  has  an  integrand 
that  has  a  one-half  power  singularity  along  variable  limits  the  variable 
limits  may  be  neglected  provided  the  finite  part  of  the  resulting 
integral  is  taken. 
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The  term  "finite  part"  is  somewhat  misleading  since  the  finite 
part  of  an  integral  can  be  infinite.  In  certain  cases  the  integral 
is  infinite  even  after  the  terms  which  approach  infinity  as  RT 
approaches  zero  have  been  neglected.. 


Scalar  Potential 


The  preceding  arguments  show  that  the  finite  parts  of  equation  (17) 
can  be  equated  to  zero;  thus. 


-2*0(x,y,z) 


n^  da  =  0 


where  the  symbol  f  before  the  integral  denotes  that  only  the  finite 
part  is  to  be  taken.  The  preceding  equation  may  be  solved  for  the  value 
of  the  potential  at  the  point  (x,y,z)j  the  result  of  this  operation  is 
given  by 

0(x,y,z)  =  2^  J  V0  -  p  .  nh  da  (l8) 

S1 

It  should  be  remembered  that  surfaces  can  exist  inside  the  forward  Mach 
cone  from  the  point  (x,y,z)  across  which  i  V0  -  0V  i  can  he  discon- 

(1  i  \  h  K 

—  V0  -  0V  —  J  .  n^  remains  continuous . 

Equation  (l8)  is  an  expression  for  the  scalar  potential  at  the 
point  (x,y,z)  in  terms  of  the  potential  and  its  derivatives  with  respect 
to  on  the  surface  S^_.  A  more  general  expression  for  the  scalar 

potential  than  that  given  by  equation  (l8)  can  be  obtained.  If  within 
the  volume  enclosed  by  the  forward  Mach  cone  from  the  point  (x,y,z)  and 
the  surface  S;l  equation  (l)  is  not  satisfied  and  V2h  operating 
upon  0  yields 


t 


V2!! L0(  =  f( 

then  equation  (10)  becomes 


(19) 


R' 


r 

V0  •  rk+ 


T 


^)da  +  fSl[^  '  ^  i)  '  ^  da 


ji  Q 


dv 


v0 


(20) 
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Provided  that  f(i,T),/;)  is  always  finite,  the  right  side  of  equation  (20) 
is  finite;  furthermore,  the  right  side  of  equation  (20)  remains  finite 
as  R'  approaches  zero.  If  in  equation  (20)  R'  is  made  to  approach 
zero  and  only  the  finite  parts  of  the  integrals  are  retained,  then  the 
resulting  expression  is 


0(x,y,z)  -  i  /  (k  Vfl  -  fa  i)  .  nh  da  -  hMlil  dv  (21) 

S1  V1 

where  vp  represents  the  volume  vQ  when  R'  is  equal  to  zero.  Equa¬ 
tion  (21)  is  equation  (58)  of  reference  4  where  p2  has  been  set  equal 
to  one.  Note  that  the  volume  integral  in  equation  (21)  has  the  appear¬ 
ance  of  the  integral  for  the  potential  resulting  from  a  volume  distribu¬ 
tion  of  sources  in  an  incompressible  flow. 

The  assumption  has  been  made  that  0  is  continuous  throughout  the 
volume  vq.  It  is  also  assumed  that  no  surfaces  exist  inside  vq 
across  which  60/Snh  is  discontinuous.  If  equation  (21)  is  applied  to 
a  volume  v,,  which  has  surfaces  across  which  0  and/or  the  derivative 
of  in  the  direction  of  nh  is  discontinuous,  these  surfaces  of  dis¬ 
continuity  can  be  removed  from  the  volume  of  integration  by  allowing  the 
arbitrary  surface  Sq  to  envelop  them  (see  fig.  3)-  For  volumes  of 
integration  where  the  surfaces  of  discontinuity  have  been  removed  in 
this  manner,  the  scalar  potential  can  be  written  as 


0(x,y,z)  = 


R 


dv 


+  sf/ 


1 

2* 


si-so 


|  v/  -  0v  i)  .  nh  da 


(22) 


where  Sq  denotes  the  surface  of  discontinuity,  and  A0  is  the  poten¬ 
tial  difference  across  the  surface  S0.  The  notation  6/8%'  is  used 
to  denote  the  operator 


2 

-P  V 


1 


s 

3T 


+  v 


6 

Srj* 


+  V 


r  C> 

3  5T 


For  the  cases  where  no  surface  of  discontinuity  exists  inside  the 
volume  v^  and  0  and  V0  are  zero  on  the  surface  -  Sq  ,  equa¬ 
tion  (22)  reduces  to 
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0(x,y,z)  =  ~  f  —  ■ dv  (23) 

V1 

From  equation  (19) 

V2h0(x,y,z)  =  f(x,y,z)  (24) 

Note  that  equation  (23)  is  a  solution  of  the  partial -differential 
equation  (24). 

For  most  problems  in  linearized  supersonic  flow,  f(i>T),£)  is  zero 
and  0  is  zero  upstream  of  the  disturbing  body.  For  such  problems,  the 
surface  Sq  -  Sq  can  be  taken  to  be  located  upstream  of  the  disturbing 
body  where  0  and  V0  are  zero.  In  this  case,  equation  (22)  reduces 
to 


0(x,y,z) 


i  a  ^  ^  a  i\ 

R  A  Sv'  ^  SrV"  R  j  da 


(25) 


If  the  surface  S0  is  confined  to  the  £  =  0  plane,  equation  (25) 
reduces  to  equation  (10)  of  reference  6.  In  this  reference  the  boundary 
conditions  for  airfoils  are  discussed. 


Components  of  Vector  Field 

Let  F  be  a  vector  which  is  finite  and  integrable  in  a  given 
volume  (denoted  by  X)  and  is  zero  outside  the  volume  X.  To  each 
point  in  the  volume  associate  the  vector 


G(x,y,z) 


(26) 


where  denotes  the  part  of  the  volume  X  enclosed  in  the  forward 

Mach  cone  from  the  point  (x,y,z). 

From  equations  (24)  and  (26),  it  follows  that  each  component  of  G 
satisfies  the  relation 


V2hGi(x,y,z) 


-2JtFi(x,y,z) 


(27) 


where  the  subscript  i  refers  to  any  component  of  the  vector  field. 
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Let  ^o(x,y,z)  be  a  scalar  and  A(x,y,z)  be  a  vector  defined  by 
the  equations 


2  =Vh  •  G  = 


F(8,n,5 )  '  Vh  p  dv 


2mA  =  V  X  G  = 


Equation  (2h)  indicates  that 


f(8,ti,5)  xVjdv 


Vh  X  (v  X  G)  =  V(Vh  •  G)  -  V^hG  (30) 

Substituting  the  expressions  for  V  x  G,  Vh  •  G,  and  V2hG  as  given 
by  equations  (2J),  (28),  and  (29),  respectively,  into  equation  (30)  and 
solving  for  F  yields 

F(x,y,z)  =  -Vf0(x,y,z)  +  Vh  x  A(x,y,z)  (31) 

Since  F  is  an  arbitrary  vector,  equation  (31)  indicates  that  any 
finite  integrable  vector  field  can  be  expressed  in  terms  of  the  gradient 
of  a  scalar  and  the  hyperbolic  curl  of  a  vector.  Equation  (31)  has  the 
appearance  of  the^  Helmholtz  theorem  (reference  3?  p.  l87)j  however, 
since  and  A  are  found  by  integration  only  in  the  forward  Mach 

cone  from  the  point  (x,y,z),  equation  (31)  hardly  seems  to  be  a  state¬ 
ment  of  the  Helmholtz  theorem  as  is  commonly  given.  The  result  given 
by  equation  (31)  was  obtained  by  Robinson  in  reference  2. 

Hyperbolic  Vector  Potential 

Equation  (31)  indicates  that  the  perturbation  velocity  vector  can 
be  divided  into  two  parts.  One  part  is  the  gradient  of  a  scalar  func¬ 
tion,  and  the  other  is  the  hyperbolic  curl  of  a  vector  function.  The 
vector  function  is  analogous  to  the  common  vector  potential  (reference  3 
pp.  10 4  and  188) $  therefore,  the  vector  function  is  referred  to  as  the  ’ 
hyperbolic  vector  potential.  Thus,  if  q'  is  the  total  perturbation 
velocity  vector,  then 


+  Vh  X  A 


(32) 
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where  0  is  the  scalar  potential  and  A  is  the  hyperbolic  vector 
potential.  The  part  of  the  velocity  vector  which  is  made  up  of  the 
hyperbolic  curl  of  the  vector  potential  is  denoted  by  q. 

By  direct  expansion  it  can  be  shown  that 


Vh  •  q'  =  Vh  •  V0  +  Vh  (Vh  X  A)  =  V2h 0  +  Vh  •  (Vh  x  A)  =  0  (33) 

Equation  (33)  indicates  that  the  hyperbolic  divergence  of  the  perturba- 
tion  velocity  vector  is  zero. 

The  vorticity  vector  is  given  by 


H  =  V  X  qT 


(34) 


Therefore,  from  equation  (32), 

H  =  V  X  (Vh  X  A) 
or 

H  =  Vh(V  •  A)  -  V2hA 

From  equations  (2d)  and  (29) *  the  divergence  of  the  hyperbolic 
vector  potential  is  zero)  thus, 


H  =  -V2!^  (35) 

Each  component  of  equation  (35)  is  a  partial -differential  equation 
of  the  form  of  equation  (24)  >  thus,  from  equation  (23)  each  component 
of  equation  (35)  has  a  solution  given  by 

,  .  i  r 

Ai(x,y,z)  =  /  - g -  dv  (36) 


where  the  subscript_  i  refers  to  any  component  of  the  vector  H.  Since 
each  component  of  A  is  given  by  equation  (36),  then 


A(x,y,z) 


1_  f 

2  «Jv1  R 


dv 


(37) 


The  velocity  vector  resulting  from  the  hyperbolic  vector  potential  is 
therefore  given  by 
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_  _  I  f  TT 

q  =  Vh  X  A  =  /  VhX^-dv 

J  V]_ 


2  fr  (y  -  *i)Hz  -  (z  -  OH, 


(39a) 


2  f  o  (z  -  £)H,  -  (x  -  !)H. 


(39b) 


B2  fr  (x  -  |)Hy  -  (y  -  T))HX 

"'-4 — — - 


where  the  subscripts  refer  to  the  components  of  the  vector  H.  The 
results  given  by  equations  (39)  were  obtained  by  Robinson  in  reference  2 


Vortex  Sheets 

If  the  vorticity  is  confined  to  a  surface  S2,  equation  (37) 
becomes 


A(x,y,z)  =  f  —  da 


Equation  (40)  is  an  expression  for  the  hyperbolic  vector  potential 
resulting  from  a  surface  of  Vorticity.  Note  that  if  the  vorticity  is 
zero  except  on  the  surface  Sg,  then  equation  (35)  reduces  to 

V2hl  =  0 

By  removing_the  surface  S2  from  the  volume  of  integration  each  com¬ 
ponent  of  A  can  be  expressed  as  (from  equation  (25)) 


■i(x,y,z)  -  j-;  f 

t>2  Xi 


•  s 

1  dnh'  R  ) 
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where  the  subscript  i  refers  to  any  component  of  the  vector  A. 
Since  each  component  of  A  is  given  by  equation  (4l),  then 

A(x,y,z)  *  5;  (i  a  -  ffl  i)  da 
Note  that  if  AA  is  zero,  equation  (42)  reduces  to 

A(*,y,z>  - 

u  S2  h 


(42) 


(^3) 


By  comparing  equations  (40)  and  (43)  it  follows  that,  on  the  surface  S2, 


H 


(44) 


Equation  (44)  indicates  that  across  a  surface  of  vorticity  the 
derivative  of  the  hyperbolic  vector  potential  in  the  direction  of  n^r 
is  discontinuous.  Thus,  a  lifting  surface  can  be  represented  by  a 
continuous  hyperbolic  vector  potential,  while  it  can  be  shown  that  a 
thickness  effect  can  be  represented  by  a  discontinuous  hyperbolic 
vector  potential.  Note  the  contrast  with  the  scalar  potential,  which 
uses  a  continuous  potential  to  represent  a  thickness  effect  and  a  dis¬ 
continuous  potential  to  represent  a  lifting  surface. 


Further  Development  of  Scalar  Potential 


The  scalar  potential  can  be  expressed  in  forms  other  than  those 
already  presented.  Equation  (8)  is  applied  to  the  region  bounded  by 
the  arbitrary  surface  Sp,  the  forward  Mach  cone  from  the  point  (x,y,z), 
and  a  second  arbitrary  surface  S3  enclosed  in  the  forward  Mach  cone 
from  the  point  (x,y,z)  and  upstream  of  the  surface-  Sp.  A  cross  section 
of  such  a  region  is  shown  in  figure  4.  The  result  of  applying  equa¬ 
tion  (8)  to  this  region  is 


1  c$T 
R  5^* 


- 1 


*  l\ 

3nj^*  R  ) 


da  + 


1  50 


,R 


-  0 


a  ±) 

5n^*  R  j 


da  =  0 


(45) 


provided  that  is  a  solution  of  equation  (l).  Note  that  the  scalar 

potential  as  given  by  equation  (l8)  is  independent  of  01  so  that  0r 
is  arbitrary  so  long  as  it  satisfies  equation  (l)  throughout  the  proper 
volume . 
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If  for  a  finite  distance  upstream  0'  is  zero  and  remains  zero 
for  greater  distances  upstream,  the  surface  S3  may  be  chosen  in  this 
region  so  that  the  integral  over  S3  in  equation  (45)  is  zero;  thus, 


s1  \R  dV 


,  6  1 
3n^*  R 


da  =  0 


Equations  (18)  and  (46)  can  he  combined  to  yield 


0(x,y,z)  = 


^-5v7  "  (0_0,) 


n  a  1  , 

T  K — rr  ^  da 


>nh*  R 


The  only  restrictions  placed  on  0T  at  this  point  are  that  it  satisfy 
equation  (l)  and  he  zero  at  a  finite  distance  upstream.  In  many  cases 

0T  may  he  chosen  so  that  or  0  -  0T  is  zero}  therefore, 

in  these  cases,  0  can  he  expressed  as 


0(x,y,z)  = 


2*JSi  R\5^  ' 


(48a) 


0(x’y’z)  =-5 i  l  »  -  f  >55-  ff  d8 


(46b) 


Equations  (48)  are  quite  useful ^  however,  remember  that  they  apply  only 
when  0*  can  he  chosen  so  that  0’  does  not  violate  any  of  its 
restrictions. 

Rote  that  equations  (48)  can  he  applied  to  problems  where  either  0 
or  80/8nh  is  given  on  the  surface  Sq.  The  application  of  these 
equations  to  most  nonplanar  problems  of  either  type,  however,  lead  to 
quite  unwieldy  integral  equations. 
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APPLICATIONS 
Planar  Problems 


Many  problems  in  linearized  supersonic  flow  deal  with  the  surface 
of  discontinuity  confined  to  a  plane  surface  parallel  to  the  x-axis. 

In  this  section  a  general  discussion  of  this  type  of  problem  is  given. 
The  coordinates  are  located  so  that  the  surface  of  discontinuity  is  in 
the  t  =  0  plane. 


The  scalar  potential  at  an  arbitrary  point  (x,y,z)  above  the 
5=0  plane  is  (from  equation  (47)) 


0(x,y,z) 


+  (0  -  0')M 


d£  dq 


(49) 


In  this  case,  the  surface  S]_  is  the  £  =  0  plane. 

If  0'(x,y,O)  is  chosen  equal  to  0(x,y, 0)  the  potential  becomes 
symmetric  with  respect  to  the  £  =  0  plane.  Thus,  for  z  =  a 

0(x,y,a)  =  0' (x,y,-a) 


and 


0z(x,y,a)  =  -0Z ' (x,y, -a) 
For  this  case,  equation  (49)  reduces  to 


0(x,y,z) 


1  r  0z(^ti>°) 

rt  1  R 

bl 


dl  dq 


(50) 


Equation  (50)  was  given  by  Puckett  in  reference  5. 

If  0'(x,y,O)  is  chosen  equal  to  -0(x,y,O)  the  potential 
becomes  antisymmetric  with  respect  to  the  £  =  0  plane.  In  this  case, 

0(x,y,a)  =  -0' (x,y, -a) 

and 

0z(x,y,a)  =  0z'(x,y,-a) 
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Thus,  equation  (49)  is  reduced  to 


0(*,y, 


a 2  f 

P  z 


dl  dt| 


Note  that  for  surface  Sp  not  confined  to  a  plane  parallel  to  the 
x-axis,  a  choice  of  0’(x,y,z)  at  the  surface  Sp  to  equal  0(x,y,z) 

at  Sp  does  not  cause  at  the  surface  Sp  to  equal 

at  Bp,  Similarly,  choosing  ^f(x,y,z)  at  the  surface 

Sp  ^equal  to  -0(x,y,z)  at  Sp  does  not  cause  x ,  y ,  z )  ap  to 

i  30(x,y,z)  .  q  nh 

equal  9J  at  Sp. 

on 


Provided  the  discontinuities  are  restricted  to  the  £  =  0  plane, 
the  scalar  potential  can  also  he  expressed  as  follows  (from  equation  (18)) 


0(x,y,z)  =-  — 


iz(  ,,  ,B2z  , 

- 5 -  +  0(  6/1,0)  V  *6  *1 


for  positive  z.  A  comparison  of  equations  (50)  and  (51)  with  equa¬ 
tion  (52)  shows  that  the  two  terms  of  the  integrand  of  equation  (52) 
contribute  equal  amounts  to  the  potential  at  any  point  (x,y,z). 

Since  the  terms  of  the  integrand  of  equation  (52)  contribute  equal 
amounts  to  the  potential  at  the  point  (x,y,z)  as  z  approaches  zero, 
equation  (52)  must  reduce  to 


0(x,y,O) 


1 

2jt 


The  preceding  equation  can  also  be  obtained  by  examining  the  limit  of 
equation  (52)  as  z  approaches  zero.  If  this  procedure  is  done  the 
entire  contribution  of  the  second  term  of  the  integrand  of  equation  (52) 
is  found  to  come  from  the  point  at  the  apex  of  the  hyperbola  formed  by 
the  intersection  of  the  Mach  cone  from  the  point  (x,y, z)  and  the 
£  =  0  plane.  Note  that  if  the  integration  is  performed  first  with 
respect  to  r\  then,  when  the  methods  of  Hadamard  are  used,  the  point 
at  the  apex  of  the  hyperbola  is  a  singular  point  and  must  be  removed 
from  the  area  of  integration  by  a  process  such  as  is  given  in  refer¬ 
ence  4,  page  147. 
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If  0Z  is  prescribed  over  the  £  =  0  plane,  then  the  potential 
is  given  uniquely  by  equation  (50).  Similarly,  if  the  potential  is 
prescribed  over  the  £  =  0  plane,  then  the  derivative  of  0  with 
respect  to  z  is  determined  over  the  z  =  0  plane.  This  result  fol¬ 
lows  from  equation  (51)  since  prescribing  0  over  the  z  =  0  plane 
determines  the  potentials  in  the  space  above  the  z  =  0  plane;  there¬ 
fore,  it  also  determines  0Z  in  the  space  above  the  z  =  0  plane  and 
the  limit  of  0Z  as  z  approaches  zero  from  the  positive  direction. 

The  question  that  arises  is  whether  0(x,y,z)  is  uniquely  deter¬ 
mined  in  the  space  above  the  z  =  0  plane  if  0  is  prescribed  over 
certain  areas  of  the  z  =  0  plane  and  0Z  is  prescribed  over  the 
remaining  areas.  If  the  assumption  is  made  that  0  is  not  determined 
uniquely,  then  at  least  two  potential  functions  satisfy  the  condition 
that  either  0  or  0Z  is  prescribed  in  all  regions  on  the  z  =  0  plane 
and  that  0  is  identically  zero  upstream  of  a  given  point.  Let  0j_ 
and  02  denote  two  potential  functions  which  satisfy  the  same  boundary 
conditions ,  and  let  0q  denote  the  potential  function  formed  by  taking 
the  difference  between  0^  and  02.  Mathematically,  the  potential 
function  0q  is  given  by 

0Q(x,y,z)  =  01(x,y,z)  -  02(x,y,z)  (53) 

Since  0^  and  02  have  the  same  values  in  certain  regions  in  the 
z  =  0  plane  then  0q  is  zero  in  these  regions.  Similarly,  since 
50Woz  and  502/oz  have  "the  same  values  in  the  remaining  regions  of 
the  z  =  0  plane,  then  c30o/az  is  zero  in  these  remaining  regions. 

The  potential  function  0q  has  the  boundary  conditions  that  either  0q 
°r  d0q/dz  is  zero  in  all  regions  of  the  z  =  0  plane  and  that  is 

identically  zero  upstream  of  a  given  point. 

Consider  the  case  where  all  the  boundaries  between  the  regions  are 
supersonic.  (The  slope  of  the  boundaries  are  such  that  the  component 
of  the  free  stream  perpendicular  to  the  boundary  is  always  supersonic.) 
The  potential  function  0q  can  be  evaluated  by  use  of  equations  (50) 
and  (51)  for  points  in  areas  which  are  far  enough  upstream  to  be  affected 
only  by  a  region  where  0  or  0Z  is  prescribed.  For  all  points  in 
these  areas  0O  is  zero  as  indicated  by  equations  (50)  or  (51).  It 
follows  from  equations  (50)  and  (51)  that  0Q  is  also  zero  inside  the 
volume  above  the  z  =  0  plane,  which  is  affected  by  these  areas  alone. 
Thus,  the  volume  where  0q  is  identically  zero  has  been  moved  downstream. 
The  same  argument  can  be  repeated  until  the  complete  z  =  0  plane  has 
been  covered.  The  preceding  arguments  cannot  be  applied  to  cases  where 
the  regions  have  subsonic  boundaries;  however,  if  it  is  permissible  to 
distort  the  boundaries  within  a  strip  of  infinitesimal  width  these  sub¬ 
sonic  edges  can  be  converted  into  supersonic  edges  by  replacing  every 
element  of  the  subsonic  boundaries  by  a  broken  line  made  up  of  supersonic 
segments.  Such  a  procedure  is  illustrated  in  figure  5.  If  the 
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assumption  is  made  that  the  subsonic  boundaries  may  be  distorted  an 
infinitesimal  amount,  then  0q  is  zero  over  the  z  =  0  plane  and  also 
in  the  space  above  the  z  =  0  plane.  Equation  (53)  now  reduces  to 

(zSpCx^y^z)  =  02(x,y,z) 


Since  and  02  are  any  'two  potential  functions  with  the  same  values 

in  certain  regions  of  the  z  =  0  plane  and  with  the  same  partial  deriva¬ 
tives  with  respect  to  z  in  the  remaining  regions,  proof  has  been  given 
that  only  one  potential  function  exists  for  which  the  potential  is 
prescribed  over  certain  areas  in  the  z  =  0  plane  and  the  partial 
derivative  with  respect  to  z  is  prescribed  over  the  remaining  areas. 

The  boundary  conditions  for  a  zero-thickness  lifting  airfoil  with 
a  given  local  angle-of -attack  distribution  are  not  of  the  type  discussed 
in  the  preceding  paragraph.  The  conditions  prescribed  in  the  z  =  0  plane 
for  this  type  of  problem  are:  The  potential  function  0  is  identically 
zero  upstream  of  the  airfoil;  0(x,y,O)  is  zero  except  on  the  plan  form 
or  in  the  wake;  the  partial  derivative  of  the  potential  with  respect 
to  z,  0Z,  is  given  on  the  plan  form;  and  0x(x,y,O)  is  zero  in  the 
wake.  The  preceding  boundary  conditions  do  not  specify  that  0  or  0Z 
be  prescribed  in  all  regions  on  the  z  =  0  plane  since  not  0  but  0X 
is  given  in  the  wake.  For  airfoils  which  have  trailing  edges  which  are 
always  supersonic,  the  requirement  that  0  be  continuous  in  the  stream 
direction  necessitates  the  potential  in  the  wake  to  have  the  value  of 
the  potential  at  the  trailing  edge  of  the  airfoil.  In  this  case,  the 
potential  function  is  uniquely  determined.  For  airfoils  which  have  sub¬ 
sonic  trailing  edges  the  Kutta-Joukowski  condition  is  generally  applied 
to  the  trailing  edges  to  determine  0  uniquely.  If  the  assumption  is 
made  that  the  trailing  edge  can  be  distorted  within  a  strip  of  infini¬ 
tesimal  width,  then  the  requirement  that  0  be  continuous  in  the  stream 
direction  can  be  used  to  determine  0  uniquely.  If  the  assumption  is 
made  that  the  subsonic  trailing  edge  is  distorted  within  the  infini¬ 
tesimal  strip  so  that  each  segment  of  each  line  element  of  the  trailing 
edge  is  always  supersonic  (see  fig.  (6)),  then  0  is  determined 
uniquely.  It  is  well-known  that  for  airfoils  with  subsonic  trailing 
edges  there  are  an  infinite  number  of  solutions  which  satisfy  the 
boundary  conditions  as  stated  at  the  beginning  of  this  paragraph.  The 
preceding  arguments  however  prove  that  there  is  but  one  solution  for  an 
airfoil  which  has  had  its  subsonic  edges  replaced  by  broken  lines  which 
are  always  supersonic.  Note  that  it  has  not  been  proved  that  the  solu¬ 
tion  obtained  by  distorting  the  subsonic  trailing  edges  corresponds  to 
the  solution  satisfying  the  Kutta-Joukowski  condition^  nor  has  it  been 
proved  that  the  solution  of  the  distorted  trailing  edge  is  independent 
of  the  manner  of  distortion. 
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Nonplanar  Problems 

The  scalar  potential  resulting  from  the  disturbances  caused  by  a 
nonplanar  body  can  be  found  from  equation  (l8)  provided  that  both  0 
and  5^/chik  are  known  on  some  surface  Si.  Unfortunately,  0  and 
6(^nh  are  not  generally  known  on  a  surface  which  fills  the  requirement 
of  the  surface  S]_;  therefore,  equation  (l8)  appears  to  have  little 
value  in  the  calculation  of  the  potential  functions  for  nonplanar  systems 
in  general.  Certain  properties  of  equation  (l8)  are,  however,  worth 
investigating. 

The  problem  of  evaluating  the  potential  on  the  upper  surfaces  of  a 
long  rectangular  body  is  discussed.  The  assumption  is  made  that  the 
body  extends  upstream  to  infinity  and  that  the  sides  are  parallel  to  the 
free-stream  direction  except  for  small  local  variations  which  cause 
small  disturbances  in  the  stream.  Figure  7(a)  shows  the  forward  Mach  cone 
from  a  point  on  the  upper  surface  of  such  a  rectangular  body.  This  fig¬ 
ure  also  shows  that  there  is  a  certain  part  of  the  surface  of  the  rec¬ 
tangular  body  in  the  forward  Mach  cone  from  the  point  (x,y,z)  that  can¬ 
not  possibly  affect  the  potential  at  the  point  (x,y,z).  If  the  surface 
S]_  in  equation  (l8)  is  taken  to  be  the  surface  of  the  rectangular  body, 
then  equation  (l8)  indicates  that  the  values  of  0  and  50/8n^  in  the 
region  which  cannot  possibly  affect  the  potential  at  the  point  (x,y,z) 
should  be  used  in  evaluating  the  potential  at  the  point  (x,y,z).  The 
only  possible  explanation  of  this  consideration  is  that  the  integral 
of  0  and  80/8nh  caused  by  the  disturbances  in  the  "blind  spot”  add 
to  zero.  This  consideration  can  be  shown  mathematically  as  follows. 

Let  denote  the  potential  function  resulting  from  the  disturbances 

inside  the  blind  spot.  From  equation  (46),  it  follows  that 

7 1  ^  %  4?  i) aa  -  °  (5M 

Equation  (54)  indicates  that  the  potential  at  the  point  (x,y,z)  can  be 
evaluated  by  applying  equation  (l8)  to  the  surface  of  the  rectangular 
body  regardless  of  the  blind  spots.  The  same  argument  holds  for  other 
bodies  with  blind  spots. 

The  preceding  arguments  can  be  clarified  by  a  simple  illustration 
of  the  affect  of  a  blind  spot.  Consider  an  infinite  rectangular  body 
such  as  shown  in  figure  7  where  the  only  disturbances  are  caused  by  a 
small  deflected  area  with  a  constant  slope  a  with  respect  to  the  free- 
stream  direction  located  on  the  lower  surface  of  the  body.  The  leading 
edge  of  the  deflected  area  is  chosen  perpendicular  to  the  free-stream 
direction  so  that  the  potential  in  the  region  not  affected  by  the 
vertical  sides  is  of  a  two-dimensional  nature.  Figure  7(b)  illustrates 
such  a  disturbing  surface. 
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by 


The  disturbance  potential  in  the  two-dimensional  region  is  given 


0(x,z) 


crv[x  -  x2  +  p(z 

- 


z2 


)] 


where  the  lower  surface  of  the  body  lies  in  the  z  =  z2  plane  and  the 
leading  edge  of  the  deflected  area  is  in  the  x  =  x2  plane. 


A  point  on  the  upper  surface  which  has  only  the  two-dimensional 
flow  in  its  forward  Mach  cone  is  illustrated  in  figure  7(c).  The  dis¬ 
turbance  potential  for  this  point  is  (from  equation  (l8))  given  by 


0(x,y,z)  =  grt 


P2ctV(|  -  x2)(z 

R3 


d?  dT] 


Upon  performing  the  indicated  integrations  the  preceding  expression 
becomes 


0(x,y,z)  =  - 
which  reduces  to 


crV 


[x  -  x2  -  P(z  -  z2)]  ffv[x  -  X2  -  p(z  -  Z2)] 


2p 


2P 


0(x,y,z)  =  0 

This  result  is  a  demonstration  that  the  disturbances  in  blind  spots  do 
.not  contribute  to  the  potential. 

The  scalar  potential  resulting  from  the  disturbance  produced  by  a 
nonplanar  body  can  also  be  obtained  by  use  of  equations  (4-7)  and  (48) 
provided  that  the  necessary  values  of  0,  50/5nh,  0'  ,  and  50'/5nh 

are  known.  In  dealing  with  planar  bodies  0'  could  be  chosen  so  that 

0  -  0’  =  0 

and  thus  equation  (47)  is  reduced  to  equation  (50).  Similarly,  0’ 
could  also  be  chosen  so  that 


50  50' 

=  3nh* 


and  thus  equation  (47)  is  reduced  to  equation  (51) .  Unfortunately,  for 
nonplanar  bodies,  choosing  0'  equal  to  0  does  not  make  50'/5nh* 
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known  as  was  the  case  for  planar  problems  and,  similarly,  choosing 
80/ Sn^*  equal  to  -50/^n^  does  not  make  01  known.  Certain  problems 
exist  in  which  01  can  be  chosen  so  that  0  can  be  written  as  a  simple 
integral. 


Intersecting  Planes 

Many  problems  concerning  nonplanar  bodies  deal  with  disturbances 
produced  by  two  intersecting  planes  parallel  to  the  free-stream  direc¬ 
tion.  In  this  section,  methods  of  solutions  for  two  planes  intersecting 
at  various  angles  are  given.  The  component  of  velocity  normal  to  the 
surface  is  assumed  to  be  known. 


Perhaps  the  simplest  case  of  two  intersecting  planes  occurs  when 
the  planes  intersect  at  right  angles.  It  is  desired  to  find  the  poten¬ 
tial  in  space  resulting  from  the  disturbances  produced  by  the  two  inter¬ 
secting  planes.  This  type  of  problem  could  represent  an  isolated  cruci¬ 
form  tail  with  supersonic  leading  edges  undergoing  various  motions. 
Problems  of  this  type  have  been  solved  in  references  11  and  12.  The 
axes  are  chosen  so  that  y  =  0  and  z  =  0  are  the  disturbing  planes 
(see  fig.  8).  When  y  and  z  are  positive,  equation  (18)  becomes 


0(x,y,z) 


1  d0(S,O+,Q 


gfy  f  r  0(s.o+.i:) 

2”  4  *3 


da  - 


1 

2jt 


I  Mi M;0+) 


da 


p Vf 

2jt  / 

JSk 


0(M,O+) 

r3 


da 


(55) 


The  surface  Sp  has  been  taken  to  be  the  disturbing  surface}  thus,  S3 
is  the  part  of  the  y  =  0  plane  (z  positive)  bounded  by  the  z  =  0  line 
in  the  y  =  0  plane  and  the  forward  Mach  cone  from  the  point  (x,y,z). 
Similarly,  S^  is  the  part  of  the  z  =  0  plane  (y  positive)  bounded  by 

the  y  =  0  line  in  the  z  =  0  plane  and  the  trace  of  the  forward  Mach 
cone  from  the  point  (x,y,z)  (see  fig.  9). 


The  assumption  is  made  that  Mi ^0+>O  and  M(  are  known 

and  that  0(t,O+,f;)  and  0(|,t),O+)  are  unknown.  The  integrals  con¬ 
taining  0(1, 0+,^)  and  0(£,r|,O+)  can  be  eliminated  by  several  appli¬ 
cations  of  equation  (46).  Equation  (46)  is  applied  to  the  volume  on 
the  left-hand  side  of  the  y  =  0  plane  enclosed  by  the  forward  Mach 
cone  from  the  point  (x,y,z),  the  y  =  0  plane,  the  z  =  0  plane,  and 
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an  arbitrary  surface  upstream  of  the  disturbance  (see  fig.  10).  The 
result  of  applying  equation  (46)  to  this  volume  is 


If  i  -  6^  V 

2 r  as  2«  J 

b5  b5 


i  f 


da  + 


If  I  MlMl^jda  +  £2  7’ 

2l7S  H  2« 


da  =  0 


The  surface  S]_  has  been  taken  to  be  the  y  =  0  plane  (z  positive)  and 
the  z  =  0  plane  (y  negative);  thus,  is  the  part  of  the 

z  =  0  plane  in  the  forward  Mach  cone  from  the  point  (x,y,z)  (see  fig.  10) 
Adding  equations  (55)  and  (56)  yields 


0(x,y,z) 


_L  /  IMLo 

2*  Jo  R  6t) 
S3  <- 


d0'(l,O-,S) 


67  j, 

2lt  ^s3  R3 


Mhikol) 

[7(?,0+,t)  -  0'(5,O-,6]da  -  ij - & - 


g f  r  Iiill 

2,1  A4  R3 


60'(^n,O+) 


pfz  f  f  0'(g,o+^) 

2,1  75  R3 


The  potential  function  0'(!,t],O  is  chosen  so  that 


U,-a,£  )  =  0(1, a, £  ) 


where  a  is  positive.  In  this  case. 
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2ir  Ja.  .3 


1 

2* 


da  - 


(58) 


Since  0'(£,'n,£)  is  related  to  0(l>'H>£)  the  only  unknown  in  the 
preceding  equation  is  0(I,T),O+).  The  region  of  integration  S5 
becomes  the  part  of  the  z  =  0  plane  (z  negative)  in  the  forward  Mach 
cone  from  the  point  (x,y,z)  obtained  by  reflecting  the  disturbing  surface 
in  the  z  =  0  plane  (z  positive)  through  the  z  =  0  plane  (see 
fig.  10) . 

The  problem  being  considered  is  one  in  which  the  normal  derivatives 
of  the  potential  function  are  known  on  two  planes  parallel  to  the  x-axis 
and  are  intersecting  at  right  angles.  The  point  (x,y,z)  has  been 
restricted  to  positive  values  of  y  and  z.  For  the  present,  consider 
the  problem  of  finding  the  potential  above  the  z  =  0  plane  when  the 
derivative  of  the  potential  function  with  respect  to  z  is  known  on 
the  z  =  0  plane,  the  derivative  of  the  potential  function  with  respect 
to  y  is  discontinuous  across  the  y  =  0  plane,  and  the  potential  is 
zero  a  finite  distance  upstream  of  a  given  point.  From  equation  (22), 
the  potential  for  this  problem  is  given  by 


0(x,y,z)  =-  — 


hMMa, . 

■K  OT) 


1  60U,q,O+) 


where  Sq  is  the  part  of  the  y  =  0  plane  above  the  z  =  0  line 
inside  the  forward  Mach  cone  from  the  point  (x,y,z),  and  S p  is  the 
part  of  the  z  =  0  plane  inside  the  forward  Mach  cone  from  the 
point  (x,y,z).  For  positive  values  of  y  and  z  equation  (59)  reduces 


NACA  TN  2641 


31 


to  equation  (58)  if  0!(|,t|,O+)  is  assumed  to  be  the  true  value 
of  0(£,*n,O+)  when  r\  is  negative.  Since  in  the  original  problem 

was  not  defined  when  T|  was  negative,  nothing  is  violated  if 
it  is  now  defined  as  being  0t(£,T|,£)  ln  the  region  where  rj  is  nega¬ 
tive  and  £  is  positive.  The  problem  in  which  the  normal  derivative 
of  the  potential  function  is  known  on  two  planes  parallel  to  the  x-axis 

and  intersecting  at  right  angles  has,  therefore,  been  changed  to  the 

problem  in  which  the  derivative  of  the  potential  has  a  known  discon¬ 
tinuity  .  across  the  y  =  0  plane  (z  positive)  and  the  normal  derivative 
of  the  potential  is  known  on  the  z  =  0  plane.  Note  that  the  poten¬ 
tial  function  still  remains  undefined  below  the  z  =  0  plane. 


Since  the  potential  function  is  undefined  below  the  z  =  0  plane, 
it  can  now  be  defined  so  that  the  resulting  potential  function  is  sym¬ 
metric  with  respect  to  the  z  =  0  plane.  Defining  the  potential  below 
the  z  =  0  plane  so  that 


0(M>a)  = 


yields  the  desired  symmetry.  The  result  of  applying  equation  (46)  to 
the  region  below  the  z  =  0  plane  inside  the  forward  Mach  cone  from 
the  point  (x,y,z)  is 


1 

2* 


1 

R 


3^z 

2 ~ 


MhlkgDda 


0(M,O+) 

R3 


da 


where  Sg  is  the  part  of  the  y  =  0  plane  (z  negative)  inside  the 
forward  Mach  cone  from  the  point  (x,y,z).  In  applying  equation  (46)  to 
the  region  below  the  z  =  0  plane,  the  surface  Sq  had  to  be  folded 
over  the  part  of  the  y  =  0  plane  (z  negative)  across  which  80/6t] 
was  discontinuous  in  order  to  be  able  to  apply  equation  (46)  to  this 
region.  Equations  (58)  and  (60)  can  be  combined  to  yield 


0(x,y,z) 


_i  r  1  50(g,o+,p  _  1  r  1  a0(6,ri,o+) 

”  *  J s3  R  ’  *  J s4  R  % 


1  r  1  a0(i,n,o+)  1  r  1S0(|,O+,£L„ 

« Ja  R  ' V -  '  *  Jq6  R  ~  da 


(61) 
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Equation  (6l)  contains  only  integrals  of  known  expressions,  and  it 
is,  therefore,  the  solution  to  the  problem  of  two  planes  intersecting 
at  right  angles  parallel  to  the  x-axis  where  the  normal  derivative  of 
the  potential  function  is  known  on  both  planes.  Figure  11  shows  the 
cross  section  of  the  distribution  of  velocity  normal  to  the  surfaces 
for  a  problem  as  represented  by  equation  (55)  and  its  solution  as  given 


by  equation  (6l).  Note  that  negative)  in  S,,  is  the 

cidfE  n  n+'l  6£  5 

reflection  of  l  in  (t)  positive)  across  the  t)  =  0  plane 

and  that  —  in  Sg  ( £  negative)  is  the  reflection  of 

~  in  s3  (£  positive)  across  the  £  =  0  plane.  This  condi¬ 
tion  suggests  that  the  result  given  by  equation  (6l)  could  also  be 
obtained  by  utilizing  the  concept  of  reflecting  surfaces. 


The  mathematical  derivation  required  for  finding  solutions  to 
problems  consisting  of  two  planes  parallel  to  the  x-axis  intersecting 
at  various  angles  can  be  reduced  by  making  use  of  the  concept  of  the 
reflecting  surfaces.  For  this  reason,  the  result  given  by  equation  (6l) 
is  obtained  by  use  of  reflecting  surfaces.  The  potential  function  can 
be  separated  into  two  parts,  0q  and  02  j  satisfying  the  following 
boundary  conditions  on  the  disturbing  surfaces: 

60j( 1,^,0+)  a0(l,T),O+) 

^  '  =  ar— 

60,(1,  o+,0 


602(M,o+) 


602(S,o+,Q  =  60(1  o+,Q 


A  cross  section  of  these  boundary  conditions  is  shown  in  figure  12. 
Only  the  potential  function  0p  is  treated  in  detail  since  the  boundary- 
conditions  for  0- L  and  02  are  of  the  same  type.  The  normal  deriva¬ 
tive  of  0i  is  zero  on  the  q  =  0  plane;  thus,  the  tj  =  0  plane  can 
be  considered  as  a  reflecting  plane.  The  potential  function  0p  is, 

therefore,  the  potential  resulting  from  a  distribution  of  60(1 ,q,0+) 

sr 

which  is  symmetric  with  respect  to  the  T)  =  0  line  and  has  the  value 
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60l(£,T|,O+) 

of  - -  when  q  is  positive.  Figure  13  illustrates  such  a 

distribution.  The  problem  of  finding  0]_  has  been  reduced  to  a  planar 
problem  which  can  be  solved  by  use  of  equation  (50). 

Equation  (50)  was  obtained  by  defining  the  potential  below  the 
z  =  0  plane  so  that  the  total  potential  function  was  symmetric  with 
respect  to  the  z  =  0  plane.  This  result  caused  the  derivative  of  the 
potential  function  with  respect  to  z  to  be  antisymmetric  with  respect 
to  the  z  =  0  plane.  Figure  14  illustrates  the  distribution  of  the 
normal  derivative  of  the  function  across  the  z  =  0  plane.  The  problem 
of  evaluating  the  potential  function  0^  ^as  "been  reduced  to  a  planar 
problem.  Similarly,  the  problem  of  evaluating  the  potential  function  02 
can  be  reduced  to  a  planar  problem.  Figure  15  illustrates  such  a  pro¬ 
cedure.  The  original  potential  function  is  the  sum  of  0^  and-  02* 
Equation  (6l)  follows  from  the  preceding  results  for  0-j_  and  02.  The 
addition  of  0^  and  02  is  illustrated  by  figure  1 6. 

The  concept  of  reflecting  surfaces  is  now  utilized  to  find  the 
potential  resulting  from  two  disturbing  surfaces  parallel  to  the  x-axis 
and  intersecting  at  an  angle  of  ^5°-  The  axis  is  chosen  so  that  the 
x-axis  lies  along  the  intersection  of  the  disturbing  surfaces  and  one 
of  the  disturbing  surfaces  lies  in  the  z  =  0  plane  (see  fig.  17). 

The  potential  function  0  is  divided  into  two  parts,  0^  and  02.  The 
boundary  conditions  on  0^  and  02  are  similar  to  the  corresponding 
potential  functions  used  for  the  disturbing  surfaces  intersecting  at  90°. 
Figure  18  illustrates  the  boundary  conditions  for  0^  and  02.  The 
surfaces  on  which  the  normal  derivative  of  0]_  is  zero  can  be  con¬ 
sidered  as  a  reflecting  surface.  This  consideration  leads  to  the  same 
distribution  of  the  normal  derivative  of  on  the  q  =  0  plane  as 

is  given  on  the  £  =  0  plane.  Figure  19  illustrates  such  a  distribu¬ 
tion.  The  problem  of  finding  0q  for  two  disturbing  surfaces  inter¬ 
secting  at  ^5°  has  been  reduced  to  a  problem  of  two  surfaces  inter¬ 
secting  at  90°.  The  solution  of  0]_  can  be  obtained  from  equation  (6l) . 
Figure  20  shows  the  surfaces  across  which  the  normal  derivative  of  0j_ 
is  discontinuous.  Since  0q  and  02  have  the  same  type  of  boundary 
conditions,  then  02  has  a  solution  as  illustrated  in  figure  21.  The 
original  potential  function  0  is  the  sum  of  0^  and  02;  therefore, 

0  can  be  found  by  considering  surfaces  of  discontinuity  as  illustrated 
in  figure  22.  The  potential  function  0  can  be  evaluated  by  use  of 
equation  (22),  because  no  surfaces  across  which  0  is  discontinuous 
exist  and  the  values  of  A  80/5nh  are  known  across  all  surfaces  of 
discontinuity. 

Another  simple  case  of  two  disturbing  surfaces  parallel  to  the 
x-axis  occurs  when  the  surfaces  intersect  at  an  angle  of  60°.  The 
potential  function  0  is  divided  into  two  parts,  and  02.  The 
boundary  conditions  on  0]_  and  02  are  similar  to  the  corresponding 
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functions  used  previously.  Figure  23  shows  a  cross  section  of  these 
boundary  conditions.  By  use  of  a  reflecting  surface,  the  function 
can  be  represented  by  the  boundary  conditions  as  shown  in  figure  24. 

The  function  0]_  is  undefined  for  240°  of  the  total  angle  around  the 
x-axis.  The  function  0p  is  defined  as  shown  in  figure  25.  Since  no 
surfaces  exist  across  which  is  discontinuous,  the  function 

can  be  evaluated  by  using  equation  (22).  Similarly,  02  can  be  defined 
as  shown  in  figure  26.  The  Siam  of  0^  and  02  illustrated  by  fig¬ 
ure  27-  The  potential  function  0  can  be  found  by  using  equation  (22). 

In  the  preceding  paragraphs,  methods  have  been  found  for  deter¬ 
mining  the  potential  resulting  from  two  plane  disturbing  surfaces 
parallel  to  the  stream  direction  intersecting  at  certain  angles.  The 
same  method  can  be  used  to  find  methods  for  determining  the  potential 
resulting  from  two  plane  disturbing  surfaces  intersecting  at  various 
other  angles. 


ROLLING  TAILS  WITH  MULTIPLE  RECTANGULAR  FINS 


The  methods  derived  in  the  preceding  section  are  used  to  find  the 
surface  velocity  potential,  the  pressure  distribution,  and  the  damping 
in  roll  of  rolling  tails  consisting  of  four,  six,  and  eight  rectangular 
fins.  For  comparison,  these  same  quantities  are  also  presented  for  the 
planar  tail  configurations  consisting  of  one  and  two  rectangular  fins. 

An  illustration  of  the  tails  treated  is  shown  in  figure  28.  The  analysis 
is  limited  to  tail  configurations  having  surfaces  of  vanishingly  small 
thickness  and  of  zero  camber.  The  investigation  is  also  limited  to  the 
range  of  Mach  numbers  for  which  the  region  of  interference  between  the 
adjacent  fins  does  not  affect  the  fin  tips. 


Tail  Consisting  of  One  Fin 

The  pressure  distribution  and  the  velocity  potential  on  the  surface 
of  rolling  tails  made  up  of  one  and  two  rectangular  fins  can  be  obtained 
from  the  results  of  reference  13.  The  pressure  and  potential  for  the 
tail  consisting  of  only  one  fin  can  be  found  by  transforming  the  axis 
of  roll  of  the  tail  consisting  of  two  rectangular  fins. 

The  tail  consisting  of  one  fin  is  divided  into  regions  as  shown  in 
figure  29(a) •  The  velocity  potential  on  the  surface  facing  the  negative 
y-directionJ is  given  by: 
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For  region  I, 


0(x^z)  =  ^ 


xz  .  .1  RM 

Tsln  i^r1- 


1  /2x  V 

3lT  - '  - b 


—  -  z 


)(z  +  f  - 1) 


For  region  II, 


0(X,z)  =  £jp 


(62a) 


(62b) 


For  region  III, 


0(x,z)  =  ^ 


xz  .  -l./Pz  l/2x 

P  Sln  V  x  +  3  VP  +  Z 


M 


(62c) 


For  region  IV  (note  that  the  potential  in  region  IV  is  the  potential  in 
region  I  plus  the  potential  in  region  III  minus  the  potential  in 
region  II ) , 


2p 


0(*,z)  =  -i 


XZ 


i(|E  -  z  -  b 
3  \P 


(TT^TTYTT)  _  ^ 

U  A  2  P/  2P 


M- 

(62d) 


The  pressure -difference  coefficient  is  given  by: 
For  region  I, 


AC  =  §p_ 

Jtpv 


(63a) 


For  region  II, 


AC 


4pz 

P  “  pv~ 


(63b) 
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For  region  III, 


8p 

"  Jt|3V 


z  sin 


-1 


For  region  IV , 


(63c) 


(63d) 


Tail  Consisting  of  Two  Fins 

The  tail  consisting  of  two  fins  (fig.  29(h))  has  the  same  potential 
and  pressure  distribution  as  a  rectangular  rolling  wing  and  can,  there¬ 
fore,  he  obtained  from  the  results  of  reference  13.  For  each  tail  con¬ 
sisting  of  two  fins  divided  into  regions  as  shown  in  figure  29(h)  the 
pressure  and  potential  in  regions  I  and  II  are  the  same  as  the  pressure 
and  potentials  in  the  corresponding  regions  for  tails  consisting  of  one 
fin. 


Tail  Consisting  of  Four  Fins 

Each  fin  of  the  tail  consisting  of  four  fins  is  divided  into 
regions  as  shown  in  figure  29(c).  The  pressure  and  potentials  in 
regions  I  and  II  are  the  same  as  the  pressure  and  potentials  in  the 
corresponding  regions  for  tails  consisting  of  either  one  or  two  fins. 

The  regions  III  and  IV  are  affected  by  the  interaction  between  adjacent 
fins.  The  potential  in  region  TV  is  made  up  of  a  combination  of  the 
potentials  of  regions  I,  II,  and  III.  Thus,  the  only  real  problem  is 
the  determination  of  the  potential  in  region  III. 

The  potential  in  region  III  is  not  affected  by  the  tip  and  is, 
therefore,  the  same  potential  as  would  be  obtained  if  the  fins  were 
infinitely  long.  With  the  coordinate  axes  chosen  as  shown  in  figure  28(c) 
the  point  (x,y,z)  is  restricted  to  values  of  y  which  are  negative 
while  the  values  of  z  are  restricted  to  positive  values.  Note  that 
for  a  tail  with  finitely  long  fins,  the  potential  at  a  point  (x,y,z) 
in  the  region  of  interaction  is  independent  of  the  disturbances  produced 
at  points  located  so  that  their  projection  on  the  yz-plane  does  not  lie 
in  the  second  quadrant.  The  general  method  previously  derived  for 
finding  the  potential  resulting  from  two  plane  surfaces  intersecting  at 
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right  angles  can  thus  he  used  to  find  the  potential  in  the  part  of  the 
region  of  interaction  which  is  not  affected  by  the  tip. 

The  velocity  component  normal  to  the  fin  in  the  z  =  0  plane  is 
given  by 


0z(x>y,O+)  =  -py 

and  the  velocity  component  normal  to  the  fin  in  the  y  =  0  plane  is 
given  by 


0y(x,O",z)  =  pz 

Figure  30  illustrates  this  type  of  normal-velocity  distribution,  and 
figure  31  illustrates  a  cross  section  of  the  surfaces  of  normal -velocity 
discontinuity,  which  previous  results  show  can  be  used  to  obtain  the 
potential  in  the  part  of  the  region  of  interaction  which  is  not 
affected  by  the  tips. 

Note  that  in  figure  32  the  discontinuity  in  the  normal  velocity 
across  the  z  =  0  plane  is  the  same  type  as  the  discontinuity  in  the 
y  =  0  plane.  Thus,  if  bhe  potential  resulting  from  this  type  of  dis¬ 
continuity  (see  fig.  32)  is  known,  then  the  potential  resulting  from 
any  combination  of  discontinuities  of  this  type  can  be  found.  The 
potential  for  this  type  of  distribution  is  denoted  by  0q.  By  the  use 
of  cylindrical  coordinates,  as  shown  in  figure  32,  the  potential  at  the 
point  (x,p,0)  can  be  expressed  in  terms  of  0q  by 

0(x,p,e)  =  0o(x,p,e)  -  0o(x,p,e  -  |)  (64) 

Equation  (64)  follows  from  figures  31  and  32. 

The  potential  function  0O  was  evaluated  by  use  of  equation  (50) 
and  is  given  by: 

For  0  <8  <-, 

—  —2’ 


0o(x,p,0)  =  - 


JtP‘ 


Xijx2  -  P2p2  +  $2p2(l  -  2  cos20)ln 


(3p 


X  +  \[x2  -  P2p2 


2x3P  cos  6  sin’^ 


(3p  cos  G 


\Jx2  -  P2p2  sin20 


2(3^p^  cos  6  sin  6  tan~^ 


x  cot  6 


x2  -  p2p2y 


(65a) 
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For  2  —  0  <  n> 


0o(x,P,0)  =  — xic2  -  P2p2  +  P2p2(l  -  2  cos20)ln  - -|3p  + 

[_  x  +  \/x2  -  P2p2 


2xpp  cos  0  sin-1 


Pp  cos  0 
:2  -  P2p2  sin20 


2p2p2  sin  0  cos  0  tan-1  /-iL-C2t— 0 


fx2  -  p2p2 


(65b) 


From  equations  (64)  and  (65),  the  potential  function  in  the  region 
of  interaction,  which  is  not  affected  by  the  tips,  is  given  by  the 


following  equation  for  —  <0  <  jt: 


0(x,p,e)  = 


-^•kp2p2(l  -  2  cos20)ln  - 


x  +  Ux2  -  p2p2 


2xPp  cos  0  sin'1 


pp  cos  0 

f2  a2  2  .  2_ 

x  -  P  p  sin  9j 


sin  6  sin-1/— gg-  .siB  0 

Wx2  -  p2p2  cos2e 


2P2p2  cos  e  sin  0 / tan-1  -f|c.ot  0|  -  tan-1  j^-gl  \  ,  (66) 

l  1A2 . 7? 


The  potential  in  region  III  of  the  fin  surface  is  a  special  case 
(0  =  of  the  preceding  equation.  Thus,  the  potential  in  region  III 
is  given  (in  Cartesian  coordinates)  by 


0(x,O",z)  =  §g£/x  sin-1  ^ 


w r  Sln  x  - pz  ln 


X  +  Vx2  -  P2Z2 
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From  equation  (67)  the  pres  sure -difference  coefficient  is  found 

to  "be 


acp  -  % f  <«» 

As  previously  stated,  the  potential  in  region  IV  (see  fig.  29)  is 
a  combination  of  the  potentials  in  regions  I,  II,  and  III.  Assume  that 
the  fins  are  infinitely  long.  In  this  case  only  two  regions  (II  and  III) 
exist  since  regions  I  and  IV  are  affected  by  the  tip.  The  effect  of  the 
tip  can  be  taken  into  account  by  adding  a  potential  which  has  zero 
normal  velocity  on  the  fin  and  the  negative  of  the  pressure  of  the 
infinite  fin  in  the  plane  of  the  fin  outboard  of  the  tip.  The  value  of 
such  a  potential  on  the  fin  is  given  by  the  difference  between  the 
potential  of  region  I  and  the  potential  of  region  II.  (This  potential 
is  only  the  effect  of  the  tip  on  a  semi-infinite  rolling  wing) .  Thus, 
the  potential  in  region  IV  is  the  potential  in  region  III  plus  the  dif¬ 
ference  between  the  potential  of  region  I  and  the  potential  of  region  II. 
Mathematically ,  the  potential  in  region  IV  is  given  by 


(69) 


From  equation  (69)  the  pres sure -difference  coefficient  is  given  by 


Tail  Consisting  of  Six  Fins 

The  pressure  and  potentials  on  the  surface  of  the  tail  consisting 
of  six  fins  can  be  obtained  in  a  manner  similar  to  that  used  for  the 
tail  consisting  of  four  fins.  The  pressure  and  potentials  in  regions  I 
and  II  are  the  same  as  the  pressure  and  potentials  in  the  corresponding 
regions  for  tails  consisting  of  one,  two,  or  four  fins.  Regions  III 
and  IV  are  affected  by  the  interaction  between  adjacent  fins.  The 
potential  in  region  IV  is  made  up  of  a  combination  of  the  potential  in 
regions  I,  II,  and  III;  therefore,  the  main  problem,  as  for  the  case  of 
four  fins,  is  the  determination  of  the  potential  in  region  III. 
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The  potential  in  region  III  is  the  same  as  the  potential  for  a 
tail  consisting  of  six  infinitely  long  fins.  The  induced  velocities 
normal  to  two  of  the  planes  of  the  fins  are  illustrated  in  figure  33- 
For  two  plane  surfaces  parallel  to  the  stream  direction  and  inter¬ 
secting  at  an  angle  of  60°,  the  potential  in  region  III  can  be  obtained 
by  a  distribution  of  discontinuities  in  velocity  as  illustrated  in  fig¬ 
ure  34.  Note  that  the  potential  in  region  III  can  be  made  up  of  a  com¬ 
bination  of  the  potentials  from  a  velocity  discontinuity  as  shown  in 
figure  35-  The  potential  from  this  type  of  discontinuity  is  denoted 
by  0q.  By  use  of  cylindrical  coordinates  as  shown  in  figure  35,  the 
potential  at  the  point  (x,p,0)  can  be  expressed  in  terms  of  by 

0(x,p,e)  =  $0(x,p,e  -  -  0o^x,p,e  -  -  0o(x,p,e  +  (71) 

Equation  (71)  follows  from  figures  34  and  35- 

The  potential  function  0q  was  evaluated  by  use  of  equation  (50) 
and  is  given  by  the  following  equation  for  0^0  *: 

$Q(x,p,e)  =  - £P-Cp-S-0 (x  -  |Bp  sin  0)  (72) 

From  equations  (71)  and  (72)  the  potential  function  in  the  region 
of  interaction,  which  is  not  affected  by  the  tips,  is  given  by  the 

following  equation  for  ^^0 

0(x,p ,0)  =  ~  2  cos^0)  +  2  cos  0  sin  0]  (73) 

The  potential  in  region  III  of  the  fin  surface  is  a  special  case 

(0  =  7^  of  equation  (73)-  The  potential  in  this  region  is  given  (in 
Cartesian  coordinates)  by 

0(x,O_,z)  =  (74) 

From  equation  (74)  the  pres sure -difference  coefficient  in  region  III 
is  found  to  be  zero. 

The  potential  in  region  IV  is  a  combination  of  the  potentials  in 
regions  I,  II,  and  III)  this  can  be  shown  in  the  same  way  as  the  poten¬ 
tial  in  region  IV  of  the  tail  consisting  of  four  fins  was  shown  to  be  a 
combination  of  potentials  from  other  regions.  Specifically,  the 
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potential  in  region  IV  for  the  tail  consisting  of  six  fins  is  the 
potential  in  region  III  plus  the  difference  between  the  potentials  of 
region  I  and  of  region  II.  Mathematically,  the  potential  in  region  IV 
is  given  by 


S2f(x,0",z)  =  ^ 


*z2/T 


l/2x 

3U 


(75) 


From  equation  (75) 


the  pres sure -difference  coefficient  is  given  by 


Tail  Consisting  of  Eight  Fins 

The  pressure  and  potential  on  the  surface  of  the  tail  consisting 
of  eight  fins  can  be  found  by  utilizing  the  potential  functions 
used  in  finding  the  pressure  and  potentials  on  the  surface  of  the  tail 
consisting  of  four  and  six  fins.  The  pressure  and  potentials  in 
regions  I  and  II  are  the  same  as  the  pressure  and  potentials  for  the 
corresponding  regions  of  the  other  tails.  The  potentials  in  regions  III, 
IV,  V,  and  VI  (see  fig.  29(e))  are  affected  by  the  interaction  between 
adjacent  fins.  Since  the  potentials  in  regions  V  and  VI  are  combinations 
of  the  potentials  in  the  remaining  regions,  the  main  problem  is  to  find 
the  potentials  in  regions  III  and  IV. 

The  potentials  in  regions  III  and  IV  are  the  same  as  the  potential 
for  a  tail  of  eight  infinitely  long  fins.  The  induced  velocity  normal 
to  two  of  the  planes  of  the  fins  is  illustrated  in  figure  3 6.  From  the 
results  for  two  plane  surfaces  intersecting  at  an  angle  of  45°,  the 
potentials  in  regions  III  and  IV  can  be  obtained  by  a  distribution  of 
discontinuities  in  velocity  as  illustrated  in  figure  37.  The  potential 
resulting  from  the  distribution  of  discontinuities  in  velocity  as 
illustrated  in  figure  37  can  be  obtained  from  a  distribution  of  dis¬ 
continuities  as  possessed  by  the  potential  function  0q  used  in  con¬ 
nection  with  the  four-finned  tail.  The  potential  function  0Q  used  in 
connection  with  the  four -finned  tail  was  evaluated  only  in  the  region 
affected  by  the  root  sections  of  the  fins.  For  the  case  of  the  eight- 
finned  tail,  interaction  occurs  between  adjacent  fins  in  regions  which 
are  not  affected  by  the  root  sections  of  the  fins.  The  potential  func¬ 
tion  0q  of  the  four-finned  tail  in  the  region  not  affected  by  the  root 
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The  potential  in  region  III  is  a  special  case 


M) 


of  equa' 


tion  (78)*  Setting  0  =  ~  in  equation  (78)  yields  (in  Cartesian 
coordinates)  . 


0(x,O",z)  =  ^(x  sin’1  ^  -  x^T s: 


-1  pz 


+  pz  tan' 


'2.XT  -  p^z‘ 


-  p^Z1 


From  equation  (79)  the  pres sure -difference  coefficient  is  found  to 
he  given  hy 


An  8pz  /  .  -l  pz 
ACP  =  w\sxri  v 


(2x2  -  P2Z2 


By  inspection  of  figures  (35)  and  (37)  the  potential  function  in 
the  region  of  fin  interaction,  which  is  not  affected  hy  the  root  sections 

of  the  fins,  is,  for  ^  <1  0  < 

0(x,p,e)  =  0o(x,p,e)  +  sz!0^x,p,e  -  -  0o(x>p>e  -  f)  -  0o{x>p’9  ~  if) 


where  0O  is  given  hy  equation  (72).  Substituting  equation  (72)  into 
equation  (8l)  yields  (remember  that  0q  is  zero  upstream  of  the  Mach 

cone  from  the  y-axis),  for  —  0  5L 


0(x,p,0) 


■  ^  ^ 


n  ,  cos  0 
m  0  +  — =-  + 

\f2j 


pp2  cos  0  sin  9  -  (1  _  2  sin20 


— cos  9  +  sin  9 
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0;  if  x  <  &=(sin  9  -  cos  0) 
x  -  ^=(sin  0  -  cos  0):  if  x  >  ^£(sin  0 

V2  ^2 


-  COS  0 


cos  0 


Oj  if  x  <  pp  sin  0 


x  -  pp  sin  0>  if  x  >  Pp  sin  0 


NACA  TN  26hl 


The  potential  in  region  IV  is  a  special  case  (o  =  of  equa- 

jt  '  ' 

tion  (82).  Setting  0  =  ^  in  equation  (82)  yields  (in  Cartesian 
coordinates) 

0(x,O",z)  =  y-[0-  -  ^2)x  +  Pz]  (83 

From  equation  (83)  the  pres sure -difference  coefficient  is  found  to 
he  given  hy 


“p  -  ^  -  VS) 


The  potential  in  region  V  is  the  potential  in  region  IV  plus  the 
difference  between  the  potential  in  region  I  and  the  potential  in 
region  II;  thus,  from  equations  (62a),  (62b),  and  (83),  the  potential 
in  region  V  is  found  to  be  given  by 


0(x,O  ,z) 


2p  nz  /  xi/^n 

*  2  (z  '  p  J 


xz  .  -1  rl 2  '  z 

+  T  sin  V — “ 
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From  equation  (85)  the  pres sure -difference  coefficient  is  found  to 
be  given  by 


8p  rtz 


m  \\[2 


z  sin"1 
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The  potential  in  region  VI  is  the  potential  in  region  III  plus  the 
difference  between  the  potential  in  region  I  and  the  potential  in 
region  II j  thus,  from  equations  (62a),  (62b),  and  (79),  the  potential 
in  region  VI  is  found  to  be  given  by 


0(x,O“,z)  = 


'^L  -x  cos-1  —  -  x\l2~  sin-1  —  ^Z 

m  x  v  \E? . 
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From  equation  (87),  the  pressure-difference  coefficient  in 
region  VI  is  found  to  be  given  by 


Discussion  of  Results  for  Rolling  Tails 


Illustrative  plots  of  the  chordwise  and  spanwise  pressure  distribu¬ 
tions  across  one  fin  for  tails  with  various  numbers  of  fins  are  shown 
in  figure  38.  Figure  39  shows  illustrative  plots  of  the  spanwise  loadings 
on  one  fin  for  tails  made  up  of  various  numbers  of  fins. 

The  potential  function  0q  used  in  finding  the  pressures  and 
potentials  for  the  tails  consisting  of  four,  six,  and  eight  fins  could 
be  used  in  finding  pressures  and  potentials  for  tails  consisting  of  any 
even  number  of  fins  provided  that  the  region  of  interaction  between 
adjacent  fins  does  not  affect  the  tip.  The  restriction  on  the  region 
of  interaction  causes  the  range  of  validity  to  decrease  as  the  number 
of  fins  is  increased.  The  range  of  validity  could  be  extended,  however, 
by  use  of  a  pressure  or* potential  cancellation  method  such  as  given  in 
references  14-  and  15. 

From  the  potential,  the  damping  in  roll  per  fin  was  calculated. 

Table  I  presents  the  results  of  these  calculations.  Figure  to  presents 
the  variation  of  the  damping  in  roll  per  fin  with  Ap  for  tails  made 
up  of  various  numbers  of  fins.  For  a  given  Mach  number  (|3  constant), 
figure  to  shows  the  variation  of  the  damping  with  aspect  ratio.  Fig¬ 
ure  bl  presents  the  variation  of  the  damping  in  roll  per  fin  with  Mach 
number  for  tails  consisting  of  various  numbers  of  fins  with  a  fin  aspect 
ratio  of  1.5. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  October  25,  1951 


46 


NACA  TN  2641 


REFERENCES 


1.  Ursell,  F.,  and  Ward,  G.  N.:  On  Some  General  Theorems  in  the 

Linearized  Theory  of  Compressible  Flow.  Quarterly  Jour.  Mech.  and 
Appl.  Math.,  vol.  Ill,  pt.  3*  Sept.  1950,  pp.  326-348. 

2.  Robinson,  A, :  On  Source  and  Vortex  Distributions  in  the  Linearized 

Theory  of  Steady  Supersonic  Flow.  Rep.  No.  9*  College  of  Aero., 
Cranfield  (British),  Oct.  1947. 

3.  Phillips,  H.  B.:  Vector  Analysis.  John  Wiley  &  Sons,  Inc.,  1933. 

4.  Hadamard,  Jacques:  Lectures  on  Cauchy's  Problem  in  .Linear  Partial 

Differential  Equations.  Yale  Univ.  Press  (New  Haven),  1923. 

5.  Puckett,  Allen  E. :  Supersonic  Wave  Drag  of  Thin  Airfoils.  Jour. 

Aero.  Sci.,  vol.  13,  no.  9,  Sept.  1946,  pp.  475-484. 

6.  Heaslet,  Max.  A.,  and  Lomax,  Harvard:  The  Use  of  Source-Sink  and 

Doublet  Distributions  Extended  to  the  Solution  of  Boundary -Value 
Problems  in  Supersonic  Flow.  NACA  Rep.  900,  1948.  (Formerly 
NACA  TN  1515.) 

7.  Martin,  J.  C.:  Retarded  Potentials  of  Supersonic  Flow.  Quarterly 

Applied  Math. ,  vol.  VIII,  no. 4,  Jan.  1951*  PP*  358-364. 

8.  Kellogg,  Oliver  Dimon:  Foundations  of  Potential  Theory.  Frederick 

Ungar  Publishing  Co.  (New  York),  1929. 

9.  Miller,  Frederic  H. :  Partial  Differential  Equations.  John  Wiley  & 

Sons,  Inc.,  1941. 

10.  Lomax,  Harvard, '  Heaslet,  Max.  A.,  and  Fuller,  Franklyn  B.:  Formulas 

for  Source,  Doublet,  and  Vortex  Distributions  in  Supersonic  Wing 
Theory.  NACA  TN  2252,  1950. 

11.  Bleviss,  Zegmund  0.:  Some  Roll  Characteristics  of  Cruciform  Delta 

Wings  at  Supersonic  Speeds.  Preprint  of  paper  presented  at 

Annual  Summer  Meeting,  Inst.  Aero.  Sci.  (Los  Angeles),  July  12,  1950- 

12.  Ribner,  Herbert  S.:  Damping  in  Roll  of  Cruciform  and  Some  Related 

Delta  Wings  at  Supersonic  Speeds.  NACA  TN  2285,  1951. 

13.  Harmon,  Sidney  M. :  Stability  Derivatives  at  Supersonic  Speeds  of 

Thin  Rectangular  Wings  with  Diagonals  ahead  of  Tip  Mach  Lines. 

NACA  Rep.  925*  1949-  (Formerly  NACA  IN  1706. ) 


NACA  TN  2641 


47 

\ 


14.  Goodman,  Theodore  R. :  The  Lift  Distribution  on  Conical  and  Non- 

conical  Flow  Regions  of  Thin  Finite  Wings  in  a  Supersonic  Stream. 
Jour.  Aero.  Sci. ,  vol.  16,  no.  6,  June  1949,  pp.  365-37^-- 

•15.  Mirels,  Harold:  Lift-Cancellation  Technique  in  Linearized 
Supersonic-Wing  Theory.  NACA  TN  2145,  1950. 


48 


NACA  TN  2641 


TABLE  I 

DAMPING- IN-ROLL  COEFFICIENT  PER  FIN 


Number  of 
fins 

3Cip  per  fin 

Valid 

for 

1 

l  +  4Ap  -  24a2P2  +  32A3p3 

24A3p3 

Ap  ^  1 

2 

1  +  8A(3  -  ^8A2p2  +  64A^p3 

48A3p3 

AP  >  | 

4 

4 

jtA3P3 

-i  +  A  +  8Ap  -  48A2p2  +  64a3P3) 

Ap  >  1 

6 

4 

A3p3 

+  — i—  (l  +  8Ap  -  48A2p2  +  64A3p3) 
9p  192  V 

ft 

8 

4 

A3P3 

—  -  \  +  —  (1  +  8a p  -  48A2p2  +  64A3p3) 

9 it  6  192  V  il 

AP  >  (2 

^NACA^ 


3.-  A  cross  section  of  the  region  in  the  forward  Mach  cone  from 
the  point  (x,y,z)  showing  the  surface  Sp  enveloping  a  surface  of 
discontinuity  S . 
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(a)  Subsonic  element  within 
infinitesimal  strip. 


(b)  Broken  line  replacing 
line  element. 


Figure  An  illustration  of  a  method  of  replacing  a  subsonic  line  ele¬ 
ment  with  a  broken  line  made  up  of  supersonic  segments. 


(a)  Wing  with  subsonic  edge.  (b)  Line  element  at  a  break  in  a 

subsonic  trailing  edge. 


(c)  Possible  distortion  of  a  subsonic  trailing  edge  at  a  break  in  a  sub¬ 
sonic  trailing  edge. 

Figure  6.-  Methods  of  distorting  a  subsonic  trailing  edge  to  determine  the 

potential  function  uniquely. 
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(b)  Bottom  view  of  rectangular  body  with  deflected  area  on  the  lower 

surface. 


(c)  Side  view  of  rectangular  body  and  the  forward  Mach  cone  from  the 

point  (x,y, z). 

Figure  7.-  Concluded. 
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(b)  Solution  to  problem  as  given  by  equation  (6l). 

11.—  A  cross  section  of  the  distribution  of  the  velocity  component 
normal  to  the  z  =  0  and  the  y  =  0  planes  represented  by  equations  (55) 
and  (61). 


6o 


NACA  TN  26Al 


Figure  12.-  A  cross  section  of  the  distribution  of  velocity  normal  to 
the  disturbing  surfaces  for  the  potential  functions  0,  0-p  and  0 


Figure  li|.-  The  normal  derivative  of  the  potential  function  across  the 
z  =  0  plane  obtained  by  applying  equation  (50)  to  a  planar  problem. 
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Figure  17.-  Position  of  coordinate  axes  for  disturbing  surfaces 

intersecting  at  45°. 


Figure  18,-  A  cross  section  of  the  velocity  distribution  on  the  disturbing 
surfaces  for  the  functions  0j_  and  02  for  the  disturbing  surfaces 
intersecting  at  45°. 
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Figure  30 A  cross  section  of  the 
normal-velocity  distribution  on 
two  plane  surfaces  representing 
the  region  of  interaction  for  a 
rolling  tail  with  four  fins. 


Figure  31.*“  A  cross  section  of  the 
velocity  discontinuity  distribu¬ 
tion  used  to  find  the  potential 
in  part  of  the  region  of  inter¬ 
action  for  a  tail  of  four  fins. 


Figure  32.-  A  cross  section  of  the  velocity  discontinuity  distribution 
associated  with  the  function  0q  for  a  tail  with  four  fins. 
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Figure  33.-  A  cross  section  of  the 
velocity  distribution  normal  to 
the  planes  of  two  fins  of  a 
rolling  tail  consisting  of  six 
fins . 


Figure  3k. -  A  cross  section  of  the 
velocity  discontinuity  distribu¬ 
tion  used  to  find  the  potential 
in  part  of  the  region  of  inter¬ 
action  for  a  tail  of  six  fins. 


Figure  35. ~  A  cross  section  of  the  velocity  discontinuity  distribution 
associated  with  the  function  0q  for  a  tail  of  six  fins. 
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Figur6  36*“"  A  cross  section  of  the  normal  velocity  induced  on  the  planes 
of  two  fins  of  a  rolling  tail  consisting  of  eight  fins. 


Figure  37.-  A  cross  section  of  the  velocity  discontinuities  which  can  be 
used  to  obtain  the  potential  in  the  region  of  interaction  between 
adjacent  fins  for  a  tail  consisting  of  eight  fins. 


Figure  38.-  Continued. 


Mach  fines 


Figure  38.-  Continued. 


Figure  38.-*  Concluded. 


Spanwise  loading  on  a  fin  of  aspect  ratio  1.5  at  a  Mach  number 
f  \f2  for  tails  consisting  of  rectangular  fins. 


Figure  1+0.  The  damping  in  roll  per  fin  for  tails  consisting  of  rectangular 
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